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~^ ABSTRACT 

We present a simple functional form for the joint distribution of i?-band luminosity and [Oil] 
ly-s I 3727 emission-line equivalent widths of galaxies, and show that this form is a good fit to the galax- 

^^ ' ies in the Las Campanas Redshift Survey. We find a relationship between [Oil] equivalent width 

^^ , W and i?-band luminosity L/j of the approximate form: {W) ~ {\QA){Lji/ Lji^^,)^^''^ , where i/?,* 

f— ^ ' is the characteristic luminosity in the Schechter function. Because this joint distribution yields 

^^ ■ information about the relationship between stellar mass in a galaxy and its recent star-formation 

rate, it can be useful for testing theories of galaxy formation. Furthermore, understanding this 
Qh' joint distribution locally will make it easier to interpret the evolution of [Oil] emission- line widths 

Q . to higher redshifts. 

^ ' Subject headings: galaxies; luminosity function, mass function — galaxies: fundamental param- 

eters — galaxies: statistics 
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1. Motivation 

Modern redshift surveys such as the Las Campanas Redshift Survey (LCRS; Shectman et al. 1996) have 
large, homogeneous sets of spectra from which one can measure star-formation indicators such as the [Oil] 
3727 A forbidden line. It is known that the luminosity function of galaxies is dependent on the emission-line 
properties of the galaxies under consideration (Lin et al. 1996; Cowie et al. 1996; Ellis, et al. 1996; Small 
et al. 1997) but the detailed relationship between these emission-line properties and galaxy luminosity has 
not been explored. Here we present a calculation of the joint distribution of i?-band luminosity and the 
equivalent width of the [Oil] 3727 line for LCRS galaxies, as well as an analytic form for this distribution 
which fits the data well. This joint distribution is a useful quantity to compare with the predictions of galaxy 
formation models {e.g., Cen & Ostriker 1998, Pearce et al. 1999, Somerville et al. 1999, Kauffmann et al. 
1999). 
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While the LCRS is the largest completed redshift survey to date, there are at least three drawbacks to 
the sample to be kept in mind. First, the survey is i?-band selected and limited by central surface brightness; 
thus, the latest type galaxies, which typically have the strongest emission lines, are preferentially excluded 
from the survey, potentially biasing our results. Second, the fits to the equivalent widths of the emission 
lines fail for some galaxies, because their spectra do not have sufficient signal-to-noise to measure the line. 
It is likely that the failure rate of the fit depends on the true equivalent width of the line, and this unknown 
incompleteness is a potential worry. On the other hand, we show below that our results are robust to the 
lower limit of equivalent widths we consider. Given the typical equivalent width errors of 2A, our results are 
most appropriate for galaxies with equivalent widths > AA, to which we limit our sample. Redshift surveys 
underway, such as the Sloan Digital Sky Survey (SDSS; York, et al. 2000) and the Two-degree Field Galaxy 
Redshift Survey (2dFGRS; CoUess 1998), will be able to overcome these two difficulties. A final problem, 
noted by Kochanck, Pahrc, & Falco (2000), is that the spectra are taken using fibers with a diameter of 
~ 3", which for typical distances of galaxies in the sample is about 4 h~^ kpc. This may cause an "aperture 
bias" which underestimates the equivalent width of emission lines at low redshift because the fiber probes 
the inner, bulge component of spirals, rather than their disks, which contain the bulk of the star-formation. 
The SDSS may be able to constrain this effect by examining the four optical colors which the survey will 
measure, and comparing the colors within fiber-sized apertures to the global colors of each galaxy. 

This paper is organized as follows. Section 2 briefly describes our method for calculating the joint 
luminosity and equivalent width function, and presents a simple fitting function based on that of Schechter 
(1976). Section 3 describes the results using LCRS i?-band luminosities and equivalent widths of [Oil] 3727 
measured by Lin et al. (1996). Section 4 suggests directions of future research. 



2. Joint Distribution of Luminosity and Equivalent Width 

We follow Sandage, Tammann, & Yahil (1979) and Efstathiou, Ellis, & Peterson (1988), maximizing 
the conditional probability that each galaxy j, given its redshift Zj, has its measured luminosity Lj and 
equivalent width Wj: 



p{Lj,Wj\zj) 



p{Lj,Wj,Zj) 



Here iminfz ) ^nd iniax(z ) &i'6 the minimum and maximum luminosities observable at redshift Zj, given the 
flux limits of the field which contains galaxy j. Wmin and T4^max sue the minimum and maximum values 
of the equivalent widths of our sample. (Lin et al. 1996 find the minimum observable equivalent width 
to be approximately constant with redshift). fgim) represents the magnitude dependence of the redshift 
completeness. The likelihood of a given model for $(L, W) is given by the product of this conditional 
probability over all galaxies in the sample. Since this conditional likelihood is independent of density, the 
normalization must be calculated separately. We use the simple estimator: 
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where V is the size of the volume probed, and (j){z) is the selection function: 
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ft is the local sampling fraction.^ Lin et al. (1996) find for the luminosity function that this estimator 
yields similar results to the minimum variance estimator of Davis & Huchra (1982) for this sample. 

We use two models to describe ^{L, W). First, we use the non-parametric form described by Efstathiou, 
Ellis, & Peterson (1988), whose extension to the two dimensional plane of L and W is trivial. Essentially, 
this method divides the (L, W) plane into bins of equal logarithmic width, and assumes the distribution 
within each bin is constant. A fast iterative method can then find the set of values which maximize the 
likelihood, and we can estimate the errors by evaluating the second derivatives of the likelihood function at 
the fitted values. 

Second, following Sandage, Tammann, & Yahil (1979), we find the maximum likelihood fit to a 
parametrized function. To do so, we parametrize the joint function as a modified Schechter function, which 
is motivated by the results below: 
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where the conditional equivalent width function is: 
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That is, at each luminosity, the equivalent widths are distributed log- normally about a mean value which 
can be expressed as a function of luminosity as: 
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(j\Y parametrizes the width of the log- normal distribution. We use this function, and maximize the likelihood 
in Equation (1) over the five parameters L*, a, Wq, iwi and A. 

For the parametric fit, we calculate the error bars using 200 Monte Carlo realizations. For each realiza- 
tion, we take the rcdshifts of all the galaxies in the actual LCRS sample to be the redshifts for the "galaxies" 
in our realization. Then, we select a luminosity and [Oil] equivalent width for each galaxy using Equation 
(4) , limiting the range of absolute luminosities for each galaxy to that which is within the flux limits at that 
redshift. Then we maximize the likelihood for this realization. This procedure allows us to examine the 
distribution of the parameters over all the realizations, and thus calculate the error bars, and to determine 
whether our method is biased. We are also able to directly compare the likelihood values for the realization 
to the likelihood value of the data sample. If the fit is consistent with the data, these likelihoods should be 
comparable; if the fit is not consistent, the likelihood value for the data will always be smaller than that for 
the realizations. 

We calculate distance moduli assuming an Einstein-de Sitter universe. We use X-corrections of the form 
K{z) = 2.51og]^Q(l -I- z) (Lin et al. 1996). Throughout, we assume Hq — 100 h km/s/Mpc with h — 1\ to 



^For a fuller explanation of the meaning of the quantities fg and ft, consult Lin et al. (1996). 



convert to other values of h, the absolute magnitude scale is shifted by 5 logj^Q h, and the luminosity function 
normalization by h^. For plotting purposes we show the luminosity function expressed per unit logarithm 
i>{L,W) = ni{ln 10) L'^{L,W). 



3. Results for the LCRS 

Here we present the joint distribution of L^ and the equivalent width of [Oil] 3727, for a sample of 
galaxies with —22.5 < Mr < —16.5 and 5,000 km/s < cz < 50,000 km/s, selected from the North and 
South 112-fiber fields in the LCRS. The equivalent widths were measured by Lin et al. (1996), by fitting 
for the position, the width, and the amplitude of a Gaussian to the continuum-subtracted spectrum near 
the predicted location of [Oil] based on the rcdshift. For about 25% of the objects, the spectra were too 
low signal-to-noise to constrain these parameters; the equivalent-width dependence of this incompleteness is 
unknown. The estimated errors in the measured equivalent widths are about 2A on average. To minimize 
the effects of incompleteness and errors, we include only measured equivalent widths > 4A in our analysis, 
leaving about 8,500 galaxies in our sample. 

Figure 1 shows the non-parametric fit as the thin solid lines with error bars. Each line shown represents 
a bin of equivalent widths, whose central value is given. Some of the lines arc offset for clarity, as described 
in the caption. Note the characteristic difference between the strong emission line galaxies, which are in 
general less luminous and have a steeper faint-end slope, and the weak emission line galaxies, which are 
brighter with a shallower faint-end slope. This result accords qualitatively with that of Lin et al. (1996) 
and those of numerous other investigations of the dependence of the luminosity function on [Oil] equivalent 
width (Cowie et al. 1996; Ellis, et al. 1996; Small et al. 1997) and on spectral type in general (Zucca et 
al 1997; Bromley et al. 1998; Folkes, et al. 1999; Loveday, Tresse, & Maddox 1999). 

We also show the modified Schechter function fit in Figure 1 as the thick solid lines for each equivalent 
width shown (again, some are offset for clarity). Apparently this model does a pretty good job, though 
it uses 6 parameters: the ordinary Schechter parameters 4><,, L*, a, plus the parameters describing the 
dependence of equivalent width on luminosity Wq^ aw, and A (which is negative, because brighter galaxies 
have smaller equivalent widths). The best- fit values of these parameters are given in Table 1. Note that it 
is approximately true from these results that 
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Also, the Schechter parameters (/>», i* and a agree generally with the results of the independent analysis of 
Lin et al. (1996), although the faint end slope here is a bit steeper. 

Table 1 also gives the errors in the modified Schechter parameters, as well as the correlation matrix 
between these parameters, determined from 200 Monte Carlo realizations, as described above. We have found 
that the bias in the maximum likelihood method is smaller than the error bars in this sample. Furthermore, 
we find that the fraction of Monte Carlo realizations which have likelihoods worse than that found for the 
data is about Pworsc ~ 0.47; this means that the likelihood for the data is comparable to the likelihoods from 
the Monte Carlo realizations, indicating that the fit is consistent with the data. 

We have experimented with performing the modified Schechter function fit with limiting equivalent 
widths between 0.5A and 9.5A, instead of the limiting value of 4A used for the results just presented. The 
parameters appear fairly robust to what this lower limit is. The largest changes are in the faint-end slope. 



which varies from a ^ —0.75 at a Hmiting equivalent width of 0.5A to a ^ —1.1 at a Hmiting equivalent 
width of 9.5A; that we measure a slightly different faint-end slope than Lin et al. (1996) is thus related to 
our choice of a limiting equivalent width of AA. Meanwhile, M, varies by about 0.15 magnitudes. However, 
the changes in the parameters which describe the distribution of equivalent widths are quite small. Wq varies 
by < 3%, aw varies from 0.85 to 0.75, and A varies from —0.45 to —0.49. This consistency simply tells 
us that the modified Schechter function is good at fitting the combination of the intrinsic equivalent width 
distribution and the incompleteness as a function of equivalent width. Nevertheless, we find it encouraging 
that the nearly same analytic form fits equally well the high equivalent width galaxies, which we are fairly 
confident of, and the low equivalent width galaxies, which may suffer from incompleteness as a function of 
equivalent width. 



4. Discussion 

We have presented a simple functional form which seems to describe well the joint distribution of 
luminosity and the equivalent width of the [Oil] 3727 emission line in the LCRS. We caution that the 
dependence of completeness on equivalent width is unknown, and further, that we have not accounted for 
the distribution of the equivalent width errors (on average about 2 A) in our analysis. Upcoming surveys 
such as the SDSS and 2dFGRS will provide larger homogeneous sets of spectra with better resolution, and 
will overcome a number of the problems encountered here. 

The joint distribution function $(Li{, W) can provide a useful tool for testing theories of galaxy for- 
mation, because the i?-band luminosity is an approximate indication of the stellar mass contained in each 
galaxy and the equivalent width of [Oil] 3727 is an approximate indication of recent star-formation in the 
galaxy. By combining hydrodynamic or semi-analytic models for galaxy formation, such as those mentioned 
above, with spectral synthesis models (Leitherer et al. 1996; Kennicutt 1998), it may be possible to place 
strong constraints on the properties of the star-formation history of galaxies. In this vein, understanding 
this joint distribution locally is also helpful in interpreting the evolution of [Oil] emission at higher redshifts 
and thus the evolution of the star- formation rate of the universe (Hogg et al. 1998). 
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Table 1. Modified Schechter Fit to LCRS Galaxies 
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Note. — Parameters of the modified Schechter function given in Equation (4) 
and correlation matrix between the parameters. The first line of the table gives the 
values of the parameters and their errors. The bottom section of the table gives the 
correlation matrix. Normalization is 0* ~ 1.34± 0.03 (xlO^^) Mpc ^^ . The errors 
and the correlation matrix were calculated using 200 Monte Carlo simulations. The 
fraction of realizations which had worse fits to the model than did the data was 
about Pworso = 0.47, indicating that the model is consistent with the data. 
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Fig. 1. — Joint distribution of luminosity and [Oil] equivalent width for the approximately 8,500 LCRS 
galaxies (in the N112 and SI 12 fields) for which we have measured equivalent widths in the range 4-100 A. 
Curves with error bars represent the results of a two-dimensional non-parametric fit based on the method of 
Efstathiou, Ellis, & Peterson (1988). They are labeled by the central value of each logarithmically spaced bin 
in equivalent width. Note the characteristic differences in M» and faint-end slope between the star-forming, 
high equivalent width galaxies, and the quiescent, low equivalent width galaxies. Smooth curves represent 
the best fit modified Schechter function of Equation (4), which appears to model the data well. Parameters 
of this fit as well as their error bars and covariances are given in Table 1. For the purposes of clarity, we 
have offset the 5.5A, 10. 5A, and 20. OA curves (for both the non-parametric and the modified Schechter fits) 
by 1.8, 1.0, and 0.3 dex, respectively. 



